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We use a numerical code to compute the density perturbations generated during an inflationary 
epoch which includes a spontaneous symmetry breaking phase transition. A sharp step in the inflaton 
potential generates k dependent oscillations in the spectrum of primordial density perturbations. 
The amplitude and extent in wavenumber of these oscillations depends on both the magnitude and 
gradient of the step in the inflaton potential. We show that observations of the cosmic microwave 
background anisotropy place strong constraints on the step parameters. 
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I. INTRODUCTION 

A period of inflation in the primordial universe pro- 
vides a causal explanation for the existence of the large 
scale structure observed during the present epoch (re- 
viewed in K|). The simplest and most natural form of 
the scalar density perturbation spectrum is the scale in- 
variant case, or Vn cx fc™ -1 with n = 1, where 1Z is the 
curvature perturbation. 

In principle, inflationary models driven by a continu- 
ously evolving scalar field have a scale dependent spectral 
index, which can be calculated using the "slow roll" ap- 
proximation |Q] . This expresses n as a function of the in- 
flaton potential and its derivatives at the instant a mode 
leaves the horizon during inflation. The inflaton evolves 
slowly, so only a small piece of the potential is "sampled" 
by the large scale structure in the present universe, en- 
suring that, if the underlying potential is smooth, n is 
not strongly scale dependent. 

Potentials with a "feature" at the value of the inflaton 
when perturbations corresponding to astrophysical scales 
in the present universe left the horizon can produce a 
primordial perturbation spectrum with significant. How- 
ever, the inflaton moves slowly and fine tuning is needed 
to put the feature in exactly the right part of the poten- 
tial. Thus, while it is possible to construct inflationary 
models with a scale dependent spectrum, these models 
are often somewhat contrived. 

However, arguing that using a feature in the inflaton 
potential to generate a complicated spectrum requires 
fine tuning assumes that the potential has just one fea- 
ture, but is otherwise smooth. Adding a large number 
of features to the potential makes it far more likely that 
a randomly chosen piece of the perturbation spectrum 
will exhibit considerable scale dependence. In particular, 
Adams et al. || showed that a class of models derived 
from supergravity theories naturally gives rise to inflaton 
potentials having a large number of sudden (downward) 
steps. Each step corresponds to a symmetry breaking 
phase transition in a field coupled to the inflaton, since 
the mass changes suddenly when each transition occurs. 



In the scenario studied by Adams et al., a spectral fea- 
ture is expected every 10-15 e-folds, so if this model had 
driven inflation it is likely that one of these features would 
be visible in the spectrum extracted from observations 
of large scale structure (LSS) and the cosmic microwave 
background (CMB). 

Motivated by the existence of models which naturally 
and generically lead to scale dependent spectra, this pa- 
per carefully examines the consequences of introducing a 
step in the inflaton potential. We focus on spectral fea- 
tures which may be observable in the large-scale structure 
or cosmic microwave background anisotropy, and there- 
fore had their origin around 50 e-folds before the end of 
inflation. 

We model the step by assuming the potential 
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for the inflaton field <j>. This potential has a step at 
4> = 0stc P with size and gradient governed by c and d re- 
spectively. For physically realistic models, inflation is not 
interrupted but the effect on the density perturbations is 
still significant. If inflation is actually interrupted the 
effect on the perturbation spectrum is severe enough to 
rule out models where this happens during the interval of 
inflation corresponding to observable scales. In order to 
accurately evaluate the spectrum, we find that we must 
evolve the evolution equations numerically, rather than 
relying on the slow roll approximation. 

Inflationary models with scale dependent spectral in- 
dices have been examined in several previous investiga- 
tions ^j-^. In particular, two recent papers, the first 
by Leach and Liddle || and the second by Leach et al. 
Q, rely, as we do, on numerical evaluations of the mode 
equation to compute the density perturbation spectrum. 
Our analysis focuses on small features in the potential. 
Conversely, Q and || examine potentials which produce 
very abrupt changes in the inflationary dynamics, includ- 
ing the temporary cessation of inflationary expansion, so 
the spectra discussed in M and |l are changed for all 
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values of k larger than some critical value. In contrast, 
the spectra we consider here are essentially unchanged 
from their form at small k once the oscillations have died 
away. Moreover, most of the models discussed in M and 
H would need to be carefully tuned in order to produce 
observable effects in the spectrum^], whereas the mecha- 
nism described by Adams et al. can alter the observable 
spectrum without fine tuning. 



II. FORMALISM 

In this section we reproduce the important equations 
governing the evolution of scalar curvature perturbations 
and gravitational waves during inflation. We use the for- 
malism developed by Stewart and Lyth |I(J] where the 
quantities of interest are the curvature perturbation 1Z 
and tensor perturbation ip. 

In the scalar case it is advantageous to define a gauge 
invariant potential 



u = —zlZ 
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where z = a<p/ H. We use the standard notation, where 
a denotes the scale factor, H the Hubble parameter, cf> 
the inflaton field, and a dot the derivative with respect 
to time t. 

The equation of motion for the Fourier components, 
u k , is 



k 2 )u k = 0, 

z 
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where the prime denotes differentiation with respect to 
conformal time and k is the modulus of the wave number 
[p~0| — |l2[| . The form of the solution depends on the relative 
sizes of k 2 and z" / z. In the limit k 2 ^> z"/z, Uk tends 
to the free field solution 
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where the normalization is determined by the quantum 
origin of the perturbations (see |l3[ ] for a more detailed 
discussion). Conversely, in the limit k 2 <C z" / z the grow- 
ing mode is 

Uk oc z (5) 
which means that the curvature perturbation, 

l^kl = K/4 (6) 



*An exception is the enhancement of the perturbations pro- 
duced just before the end of inflation, which might lead to 
the formation of primordial black holes These have ob- 
servable consequences which are distinct from observations of 
large scale structure. 



is constant in this regime. The z" / z term can be writ- 
ten as 2a 2 H 2 plus terms that are small during slow roll 
inflation, so that the first regime applies to a mode well 
inside the horizon with k 3> aH, and the second to super- 
horizon scales when k <C aH. 

The spectrum V-jz(k) is defined in the usual way as 
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and is given by 
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The mode equation for gravitational waves is 
k 2 --)v k = 0, 



(7) 
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where i>k = aipk- 111 slow roll inflation a" /a ~ 2a 2 H 2 and 
the behavior of v k is again characterized by whether the 
mode is inside or outside the horizon: 

v k -> -f== e- tkT as aH/k -► 0, (10) 

(11) 
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The power spectrum of gravitational waves V g (k) analo- 
gous to Eq. (§) is 



k 3 



(12) 



III. NUMERICAL SOLUTION 

Normally, the perturbation spectra of inflationary 
models driven by a continuously evolving, minimally cou- 
pled scalar field can be calculated using the slow roll ap- 
proximation. However, when the potential has a sharp 
feature, its derivatives with respect to <f> and the time 
derivatives of the field need not be small. Consequently, 
we evolve the full mode equation numerically, without 
any approximations other than those already implicit in 
the use of perturbation theory. 

In Eq. (^), the mode function is expressed in terms of 
conformal time. The intrinsic time-scale of the dynamics 
is not constant in conformal time, so we shift the inde- 
pendent variable to a — log a, facilitating the numerical 
integration. With this replacement, the system of equa- 
tions we are to solve is 
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where the subscript a denotes differentiation. To com- 
pute the spectrum, we repeat the integration for many 
values of k. 

In general, u has two distinct solutions since it is a 
second order linear differential equation, and we must 
choose the combination which guarantees that the mode 
equation has the limiting form, Eq. (Q). We impose the 
initial conditions when the mode is far inside the hori- 
zon assuming that the conformal time r is zero, which 
amounts to an irrelevant choice of phase. Consequently, 
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Rather than work with complex co-efhcents in the nu- 
merical code, we define two orthogonal solutions, u\ and 



ui, such that 



l k\ T=0 



du l 
da 



= 0. 



U k\ T = = - 



du 2 k 



da 



At any subsequent time Uk is thus 
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We start the evolution by evolving the two background 
equations until any initial transient solution has died 
away but the mode is still well inside the horizon. We 
then identify the two orthogonal solutions that contribute 
to Uk and extract the coefficients in (p2|). This ensures 
that an initial transient contribution to the background 
dynamics cannot contaminate the initial values of u and 
u a . Finally, to compute the spectrum, we need the 
asymptotic value of \u/z\, and we find this by continuing 
the integration until the mode is far outside the horizon 
and this value is effectively constant. 

The numerical integrations are carried out using the 
Bulirsch-Stoer algorithm jL4j, and we check our calcula- 
tions by ensuring that the results are independent of the 
distance inside the horizon where we apply the normal- 
ization, and the distance beyond the horizon where we 
evaluate the asymptotic value of \u/z\. 



IV. INFLATIONARY POTENTIAL WITH A STEP 

Figure 1 shows the power spectrum for the potential 
of Eq. (Q) with c = 0.002, or a 0.4% change in the am- 
plitude of the potential. The most striking aspect of the 
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FIG. 1. The scalar and tensor power spectrum for 
c = 0.002 and d = 0.01. The z" /z term for these parame- 
ters is shown in Figure 2. 

scalar spectrum is the scale dependent oscillations. Even 
with this small change in the amplitude of the inflaton 
potential the oscillations last for two decades of k and, 
at their peak, change the amplitude of the spectrum by 
a factor of 3. We have set the position of the step so that 
the scale where the oscillations begin, k\ ow , is probed by 
observations of the galaxy correlation function and the 
anisotropy in the cosmic microwave background. 

Before we look at the the origin of the oscillations in the 
scalar spectrum it will be helpful to have a picture of how 
inflation proceeds when there is a step in the potential. A 
general, qualitative analysis of the spectrum produced by 
a "feature" in the potential is given by Starobinsky JL5| . 
For the specific model we are considering here, we can 
understand the numerical results as follows. Energy con- 
servation requires that the change in the inflaton kinetic 
energy term cannot exceed the change in the potential 
energy so, if we are originally well inside the vacuum- 
dominated regime, a small change in the amplitude of 
the inflaton potential cannot suspend inflation. The evo- 
lution of a in Figure 2 clearly shows that the expansion is 
always accelerating. However the z" / z term, also shown 
in Figure 2, determines the growth of the scalar pertur- 
bations and is very different from 2a 2 H 2 . It first grows 
in magnitude as the inflaton field accelerates and then 
drops to a large negative amplitude as the field slows. 
However, the tensor power spectrum is unaffected since 
a" /a remains constant throughout the step. 

To understand the scalar power spectrum we begin by 
considering the evolution of a particular scalar mode. 
The evolution is governed by the competition between 
the k 2 and z" / z terms. A step in the potential of the 
magnitude we are interested in only has a lasting ef- 
fect on k modes within the horizon, and not on modes 
which are already well outside the horizon. That is the 
lowest wavenumber affected is approximately given by 
Mow ~ d-ff|step- Moreover, from the form of the k 2 — z" / ' z 
term in the mode equation, we can see that the range of 
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FIG. 2. Evolution of z"/z and a for c = 0.02 and d = 0.01 
with the number of e-folds of inflation, N. We have set N = 
at the step in the potential. 
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FIG. 3. Evolution of the independent modes u\ and u| 
(with initial conditions for it^ and u| given in Eq. (|lg|-pl| )) 
and the linear combination of their amplitude, Eq. (E2I) for 
k = 0.3. 

k affected by the step will scale roughly with the square 
root of the maximum value of z" / z in the region of the 
step. 

In Figure 3 we show the evolution of u\, uf, and uu 
for an intermediate wavenumber in the range of k af- 
fected. The rise in z" / z introduces a brief interlude of 
growing mode behavior into the oscillatory regime. The 
subsequent interval where z" / z is negative causes the 
amplitude of Uk to briefly resume its oscillatory behav- 
ior. Finally, when the inflaton field resumes slow rolling, 
the oscillations leave the horizon with altered phase and 
increased amplitude. Both of the two initially indepen- 
dent solutions are affected similarly, as they now have the 
same phase, and the amplitude of their linear combina- 
tion oscillates. In other words, the presence of the step 
introduces a boundary condition which selects a solution 
with an oscillating envelope in contrast to the uncon- 
strained plane wave solution with constant envelope seen 
at small fc. 
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FIG. 4. The effect of changing c and d on the scalar 
perturbation spectrum. 

As in the case of a featureless inflaton potential Uk 
obtains a growing mode solution once it is outside the 
horizon. However the asymptotic limit reached by the 
curvature perturbation |7£k| depends on the oscillation 
phase of the mode at horizon exit, so that |7?-k| oscillates, 
with maxima corresponding to the modes which exit at 
an extremum. The proper time interval between the step 
and when the mode with wavenumber k exits is approx- 
imately At ~ 1/afflstep - 1/aH | ex it = 1/faow - 1/k and 
in this time the amplitude of the mode will have under- 
gone l/7r(fc/fei ow — 1) oscillations. Thus the period of the 
variation in 1 72,^1 is approximately 7rfci ow . 

For higher wavenumbers the effect of the z" j z term is 
smaller, the amplitude of their oscillation is not increased 
and the two modes are not set exactly in phase with each 
other. However their phases are still altered so that the 
linear combination itk oscillates, but with a diminished 
amplitude compared to the lower k modes. 

The magnitude of z" / z depends on both the param- 
eters c and d in the potential, and in a well motivated 
model these will be determined by from particle physics. 
Alternatively, given accurate observations of the CMB 
and LSS, it may is possible cosmological constraints on 
the values of these parameters, and in the next section 
we examine the observable consequences of a scale de- 
pendent primordial spectrum. 

V. OBSERVABLE SPECTRA 

Adams et al. Q attempted to recover the primordial 
perturbation spectrum from the APM survey power spec- 
trum using the relationship between the spectrum of 
mass fluctuations today and the primordial spectrum 

Vs^VnT 2 (k)i^—j , (23) 

where T(k) is the matter transfer function that tracks 
the scale dependent rate of growth of linear perturba- 
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FIG. 5. The CMB angular power spectrum for the pri- 
mordial spectra shown in Figure 4. The normalization in 
each case is to COBE. The data are from COBE (circles), 
BOOMERanG (squares) and MAXIMA (triangles). 
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tions and depends on the dark matter content of the 
Universe. Assuming a CDM dominated Universe the 
primordial n(k) could be extracted from the three di- 
mensional PAPivi(k) inferred from the angular correla- 
tion function of galaxies in the APM survey A 
departure from scale invariance was found in the range 
k ~ (0.05 — 0.6)hMpc~ 1 . This feature has been noted 
by jlfj] and in the power spectrum of IRAS galaxies (jl7| . 
Adams et al used the n(k) they had extracted to predict 
the photon power spectrum and found that the height of 
the secondary acoustic peaks was suppressed by a factor 
of - 2. 

Recently, a number of groups have revisited this anal- 
ysis motivated by the recent Maxima and BOOMERanG 
observations which show an anomalously low second 
peak. Barriga et al performed a x 2 analysis of the COBE 
and BOOMERanG CMB data considering a simple step 
in the primordial perturbation spectrum (no oscillations). 
Griffiths et al added a Gaussian bump to the primordial 
spectrum and performed a similar exercise. Both groups 
found support for a spectral feature. 

We leave a \ 2 analysis for a forthcoming paper and 
include here for orientation the cosmic microwave back- 
ground power spectrum and matter power spectrum pre- 
dicted for the range of primordial spectra shown in Fig- 
ure 4. We use the Boltzmann code CMBFAST @ to 
calculate the CMB angular power spectrum. We use 
the less fashionable sCDM as our background cosmology 
(^cdm = 0.95, fie = 0.05, h = 0.5 ) as our motivation 
is to show the effect of the primordial density perturba- 
tion oscillations rather than find the best fit. The CMB 
angular power spectrum is shown in Figure 5 along with 
observations from COBE (the uncorrected COBE DMR 
points from |IJ), BOOMERanG §| and MAXIMA @. 
It is clear that a good fit to the data can be found, but 
that the amplitude of the step and its gradient are con- 
strained to be small by the observations. 



The matter power spectrum is shown in Figure 6 along 
with the linear power spectrum generated from the PSCz 
catalogue p4| . The theoretical spectra suffer from the 
sCDM problem of too much power on small scales how- 
ever these spectra are given for indicative purposes and 
are by no mean best-fit spectra. 

VI. CONCLUSIONS 

We use a numerical routine to accurately calculate 
the primordial density spectrum predicted by a physi- 
cally motivated inflaton potential with steps in it due 
to symmetry breaking during inflation. The step in the 
potential induces oscillations in the density perturbation 
spectrum whose magnitude and extent is dependent on 
the amplitude and gradient of the step. 

We have restricted our attention to a generic step, 
demonstrating that even a small feature in the potential 
can cause significant changes to the spectrum of large 
scale perturbations. Moreover, we have presented a de- 
tailed account of how a feature in the potential modi- 
fies the observable spectrum. In the light of our calcula- 
tions, we believe that tight cosmological constraints can 
be placed on the size of any feature in the potential, and 
thus on the particle physics model which produced it, and 
intend to return to this problem in future work. Con- 
versely, if the density perturbation spectrum extracted 
from future measurements of large scale structure and 
the cosmic microwave background turns out to be incom- 
patible with a smooth initial spectrum, the mechanism 
proposed by Adams et al. provides a natural mechanism 
for injecting significant scale dependence within the con- 
text of inflation. || 
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